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We present a protocol to generate and control quantum entanglement between the states of 
two subsystems (the system S) by making measurements on a third subsystem (the monitor M), 
interacting with S. For the sake of comparison we consider first an ideal, or instantaneous projective 
measurement, as postulated by von Neumann. Then we compare it with the more realistic or 
generalized measurement procedure based on photocounting on M. Further we consider that the 
interaction term (between S and A4) contains a quantum nondemolition variable of S and discuss 

■ the possibility and limitations for reconstructing the initial state of S from information acquired by 

' photocounting on Ai. 
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I. INTRODUCTION 



Information processing can be largely improved when quantum properties are used for encoding both, bits and 
channels jlj] . While bits are encoded in mutually orthogonal states of a quantum system, quantum channels use the 
ability to set systems in entangled states. Entanglement of states is a characteristic quantum correlation that, in 
principle, can be produced in post interacting quantum systems @. To use this quantum resource for information 
processing one has to be able first to produce and then to control the amount of entanglement of a finite number of 
quantum systems. A fundamental open problem in quantum information is the characterization and classification of 
mixed entangled states of a multipartite systems || . Nowadays the most accessible and controllable source of quantum 
entanglement has been the electromagnetic field, through parametric down-conversion processes in non-linear crystals 
J|,[|. Recently, internal atomic states entanglement have also been considered in distinct experiments 

We propose a consistent scheme for generating and controlling entangled states of two subsystems (that we call 
S): (i) If the subsystems do interact, their initial states should be adequately prepared such that the interaction Vs 
does not entangle their states during the time evolution, (ii) A third quantum subsystem, the monitor A4, is coupled 
(through Vsm) to S; M should be the only subsystem responsible for entangling the states of S, thus formally, 
[KssKsm] = is a necessary condition. Then follows an operational procedure or protocol: (iii) First one choose 
an observable K. of M. and after an elapsed time t, from the beginning of the interaction (between S and Ai), it is 
measured and the interaction is turned off; the eigenvalue outcome determines the entangled state in which S is left. 

In an ideal projective measurement the measured eigenvalues k occur with a certain probability, so the entangled 
state of S obtained through reduction, cannot, in principle, be chosen. However, if the experimentalist is able to 
control the outcomes of A4 to be read, then, controlling the degree of entanglement in S becomes possible and 
the protocol becomes feasible. Interestingly, in a realistic measurement, when quanta are counted, it is possible to 
control the outcome of the monitor: One turns on the interaction between S and M. and when k quanta are counted 
then one knows in which state the system S is left. However, the necessary interval of time for counting k quanta 
is probabilistic: if the experiment is repeated, the same number of quanta may be registered within another time 
interval. Thus, for reproducing the same state one should repeat the experiment such that the counted photons be 
the same within the same time interval. The present status for generating entangled states is very different from this 
protocol, the experiments are based on projective measurements [0J^]. 

The realization of the proposed scheme and protocol are considered in the following physical system: the subsystems 
constituting S may consist of two interacting (but not necessarily) electromagnetic (EM) fields, modes A and B, 
coupled to the monitor A4, a third EM field, the mode C. The entanglement in S is created by a continuous destructive 
photocount on C and the control is fulfilled by turning off the interaction when a certain predetermined number k 
of photons become registered, thus the system S is left in an entangled state which is essentially characterized by k. 
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This proposal is detailed in the following sections: In Sec. II we present our model and write the Hamiltonian for the 
system 5, the monitor M. and their interaction, we also give the time-dependent state vector of the whole system. 
In Sec. Ill we describe the measurement process in two different approaches for sake of comparison, the ideal and 
the realistic: (i) The ideal or instantaneous projective measurement on field C, assumes the statevector reduction by 
projection as postulated by von Neumann, leading to an entangled (regarding the fields A and B) pure states, (ii) 
Then, more realistically, we consider the measurement as a sequential photocount process, on field C, as proposed 
by Srinavas and Davies JToj] . This leads to an entangled density operator for fields A and B. We show that only the 
second approach allows controlling the degree of entanglement of the AB fields, depending on the detector counting 
rate (7), the number of counted photons (k) and counting time. In Sec. IV we show that due to the nondemolishing 
character of the coupling between S and M. , the total photon number N of S can be inferred (without altering N) by 
averaging over counted photons of mode C. We also analyze the information one gets about the initial state of S from 
the counting process on M.\ examples are presented and discussed. In Sec. 5 we present a summary and conclusions. 



II. MODEL 



The problem of production and control of state entanglement between two subsystems is based on the physical 
paradigm of four interacting EM fields. Couplings of EM modes are made possible in nonlinear media and phenom- 
ena such as parametric down and up conversion appearing when a response of second order nonlinearity in crystal 
polarization is present and four- wave mixing (third order nonlinearity) occurs in a Kerr medium fll|| . Also, recently 
a strong field-field interaction of few photons, induced by non-resonant interactions between fields and a Cs atom, 
was observed experimentally [ jl2| . This observation led to a proposal for attaining high-nonlinearities with single 
atoms fl3|| . The dynamics of the fields here considered consists of two processes: (i) a second order nonlinear process 
coupling the modes A and B (5), assisted by a classical pump field 0Jl^l and (ii) a four- wave mixing, with modes, 
A, B and C (treated as quantized fields) coupled to a fourth classical intense field. The system dynamics is described 
by the Hamiltonian 

H = hw a a)a + hiubtfb + huj^c + HX {a%e lut + atf er lvt ) 

+h X {<Ja + tfb)(ce- w,t + c s e iv ' t ) . (1) 

The total number of photons of modes A and B, N = h a + fib — a^a + b^b is a quantum nondemolition (QND) variable. 
This is a quite important feature because while measuring (destructively) the mode C, an inference can be made on 
N 2 , without loosing or altering a single quanta of modes A and B. It has been shown that the coupling between 
A and B as in Hamiltonian ([[]) displays several interesting features (i) It leads to a complete states swapping 

(information exchange) even at constant mean energy, (ii) If the states are initially not entangled, the interaction 
will produce entanglement only if one of the modes is prepared in a nonclassical state; otherwise, if both modes are 
initially prepared as a direct product of coherent states the interaction will not change this character in the course of 
their evolution J|. 

In the interaction picture, Hamiltonian (Q) is written as 

Hi = h\ (a f 6 + atf) + h x {a)a + tfb){c + c f ) , (2) 

where, resonance conditions, v = uj a — lu^ and v' = uj c , must be satisfied in order to eliminate the explicit time 
dependence present in ([!]). 



Next, suppose the monitor mode C (the meter in the terminology of |17|) is prepared in the vacuum state, so, 
during its evolution it could absorb energy only from the intense classical field. The fields A and B are considered, 
for the moment, to be in an arbitrary state, that we write as an expansion in the number states basis \m) a and \n) b, 
\iPab{0)) = Tl m n C m ,n\m, n) (from here on we omit the subscript AB). At time t the evolved state of the AB system 
is represented by 

!*(£)) =e- 4A *( Q+b+Qbt )^C m ,„|m,n)(g)e- 4Xt(m+ ™ )(c+ct) |0) c . (3) 

m,n 

Noting that e -*xt("i+™)( c + c ' ) |o) c = | — ixt(m + n))c is a coherent state restricted to the imaginary axis, state ||) can 
be written as 

!*(*)> = c ™,nUt\m, n>® I - ixt(m + n)) , (4) 

(here on we omit the subscript C) where Ut = e ~ lXt ( a b+ab ). I n the next section we discuss two forms of measurements 
on mode C. 
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III. MONITORING AB FIELDS TROUGH MEASUREMENT ON C 



A. Instantaneous projective measurement 



An ideal or instantaneous projective measurement (PM) of a system £ is associated to a chosen observable A (or a 
set of commuting observables) of this system. If at time t one of the eigenvalues A m of A is realized, then the system 
state |^(t)) is reduced to its corresponding eigenvector |m). The probability for that realization is |(\Ij(i)|m)| 2 . This 
is known as a measurement of the first kind, as postulated by von Neumann. In our model, if at time t it is found that 
field C contains exactly k photons, the state is automatically projected on to the number state \k). As a consequence, 
the AB joint state 'reduces' instantaneously to the new state 

(fe) (t) _ (k\p(t)\k) 
PabW ~1r AB [(k\p(t)\k)] (5) 

where p(t) = \^ (t)) (t)\ and Tr ab stands for the trace operation on the A and B fields operators. The probability 
that at time t the field C has k quanta is P(k,t) — Tr ab [(k\p(t)\k)]. For the state given by equation (Q) this 
probability is 

P(M) = [ ^ r ^AB{0)\UlN 2k e-^ 2fl2 U t \^A B m 

= [ ^{*AB(0)\N 2k e^ 2 * 2 \*AB(0))= £ \C m , n \ 2 ±( X t(m + n)) 2k e-^ m+ ^ 2 , (6) 

m,n— 

the second equality follows because [Ut,N] = 0, so the counts are independent of the evolution of AB modes, they are 
not affected by the evolution of the AB system. The evolved state is reduced to 



pf B (t)=P-\k,t)±( X t) 2k N k e-^ 2 * 2 U t \* (7) 
or in terms of the statevector (excepting a phase factor), 



(t)) = = N k e-iW* 2 U t \* AB (0)) , (8) 

the purity of the state is maintained, however, entanglement is created even if initially the state I^abW) is factorized. 
The mean photon number of field C is closely related to the mean (N 2 ) of AB fields, 

OO oc 

k = J2kP(k,t) = ( X tf £ \C m , n \ 2 (m + n) 2 = ( X t) 2 {N 2 ), (9) 

k— m,n— 

so, a measurement on C allows to infer the mean squared QND variable with a proportionality factor that goes as t 2 . 
The variances are related as 

Var(fc) -k= ( X i) 4 Var(A 2 ) , (10) 

where Var(x) = (x 2 ) — (x) 2 . 

If fields A and B are initially prepared in number states, |^ab(0)) = \m,n), the reduced state of the AB field 
becomes independent of k, evolving freely as 

\*AB(t)) = U t \m,n) , (11) 

thus, not feeling at all the presence of field C; so, any entanglement will only arise from the interaction between A 
and B and not from a measurement on C. The probability for measuring k photons in C at time t will depend on m 
and n as (m + n) 2 , the probability distribution being Poissonian 

P(k,t) = A \( X t) 2 (m + n) 2 ] k e -M 2 (™+«) 2 ] . (12) 
fc! 
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The mean photon number of C Eq. (||) gives 

k = ( X t) 2 (N 2 ) = ( X t) 2 (m + n) 2 , (13) 

whereas 

Var(fc) -k = (xt) 4 Var(7V 2 ) = . (14) 

Inversely, if one ignores in which states the AB modes were prepared, but however verifies, through a measurement 
on C, that Var(fc) = k, then one immediately infers that the AB system was prepared in some number state \m,n). 
If the modes A and B are prepared in eigenstates of the QND variable the subsystem AB evolves independently of C. 

Now, if the initial states of the subsystems A and B are prepared as a direct product of coherent states, \^ab(0)) = 
a) ® = \a, (3), the unitary evolution does not entangle the states of A and B, each one continues its evolution as 
such, 

|*W (i)> = ^^ N k e-^ 2 * 2 \a(t),f3(t)) (15) 

where Ut\a,/3) = \a(t), (3(t)), with a(t) = a cos Xt — i/3 sin Xt and (3(t) = (3 cos Xt — ia sin Xt (note that |a(i)| 2 -|-|/3(i)| 2 = 
|a| 2 + |/3| 2 = Fisa constant of the motion). Although the measurement on C affects the dynamics of the AB system 
by entangling the states, this very ideal measurement does not allow controlling the entanglement dynamics of the 
AB system, since the outcome of the eigenvalue k is probabilistic, according to the distribution (|^). 
As in the previous case, the mean photon number of C will depend on time as t 2 , 

k = ( X t) 2 F(F+l) , (16) 

the variance is given by 

Var(fc) = ( X t) 4 (AF 3 + 6.F 2 + F) . (17) 
and F can be obtained from the measurements on M., 

F = [(Var(fc)/( X t) 4 - 2k/( X t) 2 ) / (Ak/( X t) 2 - l)] . (18) 
where the right hand side, calculated from the experiment should be time- independent. 



B. Generalized measurement by photocounting 

Experiments involving counting are not instantaneous and far from the von Neumann idealization. For a realistic 
counting measurement, one should consider that (i) For a given time interval t the counting of k quanta occurs with 
some probability, or, it is not likely that exactly k quanta are counted in a predefined time i, actually, there will be 
a distribution of time intervals, (ii) More importantly, one has to consider that when one photon is counted the EM 
field will have one photon less. The dynamics of this kind of process has to be treated as a dissipative continuous 
measurement; this subject was well addressed by Srinavas and Davies [|l0) and applied to several situations in [|l7]-[l9). 
We follow closely these references, describing the continuous photocount measurement in the formalism of operations 
and effects Q. 

The count of k photons from the monitor mode in a time t is characterized by the linear operation N t (k), acting 
on the state of the system, 

where p(0) is the state of the ABC system prior turning on the counting process and P(k,t) = Tr[N t (k)p] is the 
probability of counting k photons in t. The linear operator N t (k) is written in terms of two other operators, S t and 

J ! 

N t (k) = f dt k f dt k -! ■ ■ ■ f 2 dtiSt-ttJStt-t^ ■ ■ ■ JS tl , (20) 
Jo Jo Jo 
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where St = N t (0) is a superoperator defined in terms of ordinary Hilbert space operators B t as 

S t p = B tP B\. (21) 
B t is a semigroup element given in terms of the generator Y as B t = e Yt . As defined in |l0) the generator is 

Y=~H-R/2, (22) 

where H is the system Hamiltonian and R is the rate operator given by 

i? = 7C t c, (23) 

and the parameter 7 stands for the detector counting rate. The theory becomes self contained by choosing J as 

Jp = -ycpc 1 * , (24) 

standing for the change of the field C due to loss of one counted photon and S t is responsible for the state evolution 
between counts. 

In the interaction picture Y becomes 

Y = -iX (a f 6 + a6 f ) - i X {^a + tfb)(c + c f ) - I c +c , (25) 

the first term contributes only to the free evolution of the AB fields as a unitary evolution of the initial state. The 
monitor field stands for the counting process, being present in the other terms, thus the linear superoperator Nt(k), 
acting on the initial state Yl m n Cm,nl TO 7 n ) ® can be expressed as 

N t (k) =U f dt k f tftfc_i • • • / 2 dtiSt-t^JStH-tu-! ■ ■ ■ JS tl , (26) 



where 



with 



Sf = e ™ • e^* (27) 

Yi = -ix(m + n)(c + J) - |c f c (28) 
% = -ix(m' + n')(c+J) - |c + c. (29) 
We have defined lit as the superoperator for the coherent evolution of modes A and B 

K t% — e -iAt(a t 6+a6 t ) # giA^at b+ab+ ) ^ /jjqs 

After doing some algebraic manipulations we find that for k counts on C the state for the ABC system becomes 
P ik) ( t ) = pK J 29 ®^ E C m , n C* m , in ,(m + n) k (m' + n') k eMA(t)+A'(t) 

' m.n 7 m' ,n ; 

+i|z(t)| 2 + ±|z'(t)| a ]« t [Kn)<mV|] [|^(t)> (^(t) |] , (31) 

where g(t) = (2x 2/ 7 2 ) (-3 + 7< + 4e~ 7 */ 2 - e"^), z(f) = {-2ixh){m + n)(l - e~ 7t / 2 ) is the label of the coherent 
state and A(t) = (— 2\ 2 h 2 ){m + n) 2 Wt — 2(1 — e 7 */ 2 )] . z'(t) and A'(t) stand for the same expressions but with 
ml + n' instead of m + n. The normalization function P(k, t) is the probability for k counted photons in time t, 

P(k,t) = [ ^t^ AB (0)\N 2k e- 2 ^ 2 \^f AB (0)), (32) 
k\ 

the same as in the ideal PM, Eq. (^), however the factor (%t) 2 is now replaced by 2g(t). An important difference 
arises here, for yt ^> 1 we have 2g(t) ~ (2x/^) 2 'jt, the time appearing linearly while in the ideal PM it is quadratic. 
The expressions for the mean and variance of k are as (|6|) and ( |l7| ) however with (2xh) 2 ~ft substituting {xt) 2 ■ 
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The mixed state p^ (t) is a post-measurement state for counting k photons. The pre-measurement state is obtained 
by summing over all the possible outcomes, ^ fc N t (k)p(0) fll9f . 

The state of the AB system is obtained by tracing over the monitor mode, 

„(*) (f) - 1 [Mt)f <rk -h(t)N' 



i2Ci"i_«.— r*/ 2 i 2 r 
e i 



-(1-e-/ ) ^[^(^^(Q))^^)!)]^-^ 2 ^, (33) 

where h(t) = (2 X 2 h 2 )["ft - 2(1 - e~T*/ 2 )], and 



{l-e—< t/2 ) 2 N*N A _ 



[42^(i _ e -7*/a)2l 
6 = — » } -N 1 6n 1 . (34) 



2=0 



Preparing the modes A and B in number states the mode C has its dynamical evolution decoupled from the AB 
system. As in the ideal PM, Eq. (0), the AB modes states evolve freely, not being affected at all by the counting 
process (due to the QND variable). 

For the modes A and B prepared in coherent states, I^abCO)) = \a, (3), the unitary evolution does not change this 
character, each state continues to evolve as such (\a(t), f3(t)){a(t), (3(t)\), however, the photocount on C affects the 
dynamics introducing a another effect: besides the entanglement, as in the ideal PM, the factor in brackets in ( |3^ ) 
mixes the states. The number of counted photons in mode C determines the selection of a specific entangled state 
of the AB field and the degree of its entanglement. This will depend on x, the counting rate of the detector 7, the 
number of counts k, the time t, and there will also be a dependence on the initial state through the total number of 
photons operator N. So, controlling these quantities entails the correlated state (|33"|). It is worth stressing that time 
interval for counting exactly k photons is probabilistic and P(k,t) is its (non-normalized) distribution function. 



C. Degree of entanglement and mixing 

It is well known that a precise measure of the degree of entanglement is not available for continuous variables 
mixed states || . We saw that a straightforward application of the photocount measurement process, may be used for 
producing and determining the degree of a fc-entangled state of AB modes. Formally, for the AB system, prepared 
initially in coherent states, the dissipative (nonunitary) character of the evolution is induced by the counting process, 
being responsible for the interplay between the entanglement, mixing and decoherence in ([33]). 

For the lowest values of k the probability distributions P(k,t) are depicted in Fig. 1 for \a\ 2 = \fi\ 2 = 5. One 
perceives that for k — the highest values of the probability occur for -C 1, however for k ^ the probabilities 
attain maximum values at different times t rn . For jt 1 one has h(t) w (x^) 2 /2, 4(x/7) 2 (l — e -7 */ 2 ) 2 w (xt) 2 an d 
2fif(*) ~ (x/7) 2 (7^) 3 /3, so state ( |33| ) does not mix substantially for small time intervals, it can be written as 

*P- 1/2 (k,t)N k \a(t),0(t)), (35) 

(here, P(k,t) — (a(t), [3{t)\N 2k \a(t), f3(t))) and, excepting for k = (P(fc,0) = Sk,o), the counting process generates 
entanglement. 

Comparing state (^) with (|l^) for (x/l)lt *C 1, we see that they are very likely, thus for small times they cannot be 
distinguished and the entanglement will be stronger the higher the number of counted photons. We give a quantitative 
picture of entanglement by calculating the excess entropy [^lj , defined as 

I = S A + S B - Sab, (36) 

where Sa (Sb) is the entropy of the mode A (B) (associated to p\(f) = TtbPab(^)) an< ^ is the entropy of the 
AB system, where 5 = 1 — Trp 2 . The excess entropy measures the information contained in the correlation between 
modes A and B. Its lower and upper bounds are obtained by the Araki-Lieb p2| inequality 

\Sa — Sb\ < Sab < Sa + Sb, (37) 

namely, 



G 



< I < 2 mm(S A ,S B ). (38) 

Remark that / measures the correlation between A and B, without resolving between classical or quantum correlation 
(entanglement). However, when the joint system is in a pure state Sab = and Sa — Sb, the inequality j38| ) reduces 
to 

< I < 2S A - (39) 

and any correlation given by / > will be due to entanglement. Thus, the measure of / together with the system 
degree of purity allows one to distinguish between classes of entangled states, even though it does not give a precise 
borderline for separability of mixed entangled states. 

The excess entropy for the system here considered is plotted in Fig. 2 as function of k for most probable time of a 
/c-event (crosses for (|33|)) and for initial times (7^ ^ 1) (filled squares for (|35|)). It is verified that the counts correlate 
the subsystems more intensely the higher is the number of counted photons. Since state (|35| ) is pure (Sab = 0), this 
correlation characterizes a maximal entanglement. When detections occur at initial time the state of the joint system 
AB is left in a pure entangled state and the degree of entanglement is directly proportional to the number of counted 
photons. The excess entropy for state ( [33] ) is calculated at times tm (for which P(k, tm ) is maximum - as in Fig. 
1) and one sees that the correlation is stronger than for state (^5|). However, as state (ra) is mixed one also need 

to measure its purity (Sab), as depicted in Fig. 3, for several values of k and at times tm ■ The combined results 
of Figs. 2 and 3 for the mixed state entanglement, show that the higher the number of counted photons, the more 
correlated becomes the AB state, and for k ^> 1 one has J« 1. However, the states more likely to occur at times 
tm^ are not pure anymore and for high k's, characterized by I < 1 and Sab ^ 1, the state (^3|) becomes classically 
correlated (separable), while for < Sab < 1 it is non-separable, which characterizes non-maximal entanglement. 

How this process happens can be represented by the following example, which exactly matches our results for both 
the excess entropy and purity. Considering the condition of separability, 

PAB = P^A ®°B> J2 Pi = 1 > ( 4 °) 

1 i 

with pab acting on Ha ^H-b] dim 7-^ = Na and dimTis = ^b-> then, 

1=1- Y,PiP) [toA'M + T±b*$<t$ - (^A^af) (Tr B a$a^)} . (41) 



^,3 



When 



one obtains 



°X = Nl and °b=n- b (42) 



(N^ 1 +N B X -(N A N B ) x ) and S AB = l-(N A N B ) 1 (43) 



1=1- 

For Na,N b — > oo (as it is for coherent states Vboth quantities go to 1. So, for a k sufficiently large the state (|33 
becomes separable. State ( p0| ) (together with (f42|)) corresponds to an equiprobable ensemble of states. 

In that way, we saw that a full range of states may be generated, with different degree of entanglement, from 
maximally entangled to separable states. So, the protocol for producing a specific state ([53]) consists in turning off 
the interaction when an a priori selected value k is attained at a time t. 

Besides, the photocount measurement theory and the specific Hamiltonian (|^) also allow to extract information 
about the AB system by counting the photons of C. Due to the nature of the coupling between the monitor and 
the AB system, the counting on C gives the maximal information about any function of TV 2 . Then it is immediate 
to put the question: How much information may the photocount distribution give about the initial quantum states 
|^ab(0))? This is discussed bellow. 

IV. PROBING LIGHT WITH LIGHT 

Here we look at the counting problem as a spectroscopic measure of the state of the AB system. Light state 
spectroscopy, or light probing via atomic deflection was discussed by M. Freyberger and A. M. Herkommer p3[. 
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Probability distribution of transversal momenta of two-level atoms deflected by an EM field in a cavity allows the 
complete knowledge of the field quantum state. Is it possible to use a similar strategy with fields only, but now instead 
of the beam deflection, the photocount playing the role revealing the AB state? 
To make this point clear let us consider the average counted photon number 



k = u{t)(N 2 ) =u(t)Y,\C m ,n\ 2 {m + n) 2 



(44) 



where u(t) = {2x/lf (-e~ 7t + 4e~ 7 * /2 + 7* - 3). For it > 1, u(t) = {2x/lf it and we can write 



k = 



(2 X /l) 2 lt\j2\C m ,n\ 2 (m + nf 



(45) 



and (N 2 ) can be inferred directly by the mean counted photons, this being more precise the larger yt is because one 
gets a linear relation between time t and k. This point was actually pointed by Milburn and Walls |17| for a single 
mode coupled to C. Note that differently from the projective measurement where k increases with i , Eq. (|l3|), in 
the non-ideal measurement k increases linearly with time. We can build up higher-order moments of k in a similar 
fashion, obtaining, 



2 -u(t)(m+ny 



du 



a u(t)(m+n) 



(46) 



or keeping only the higher power in 



{2 X hf {it) 



^>> 1 for conveniently chosen time t, we have 



k r = 



(2 X /lf (it)} r J2\C m ,n\ 2 (m + n) 2r . 



(47) 



Calling 



K = k/ (2 X /lf(lt) 



(48) 



multiplying both sides of the above equation by (— l) r x 2r /r\ (0 < x < 2ir) and summing over r we get a new function, 
which a Fourier expansion of the squared moduli coefficients, 



°° (-lYx 27 "^ 

H(x) = ^2 j = X! |C™.™| 2 cosir(m + n), 



(49) 



which is always bounded, with values in the interval [—1,1] since J2 m n I Cm,™ 1 2 = 1. The left hand side is calculated 
from the experimental photocounts, thus a specific relation between a particular set of the |C TO! „| may be inferred 
in certain cases, as to be see below. Doing an inverse Fourier transform we get 



C{j) = - / cos jxH(x)dx = Y |C mj _ m | 2 . (50) 

which is the only information on the initial state of the AB system one gets by photocounting on C. If the fields A 
and B arc initially disentangled, |C mi „| 2 = |^4 m | 2 |i? n | 2 , then in order to determine the coefficients of, for instance, the 
initial state of mode A, |^4 m | 2 , one should set the initial state of B in the vacuum state (or any other number state), 
since |i?o| 2 = 1 and all other coefficient being zero, so one obtains 

\A m \ 2 =C(m). (51) 

This strategy has its limitation since only the moduli of the coefficients can be obtained, being the maximum infor- 
mation about the initial state of field A that one can obtain by counting on C. 
Particularly interesting states to be addressed are the initially entangled states. 

(a) We first consider the case when both modes are prepared in the superposition of perfect anti-correlation 
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JV 



\*Ab) = C N-n,n\ N ~ n > n ) 



n=0 



defined as a limited sum of the photon number state, and X^=o 1^ 
of state is for entangled qubits (dim[7i^ ® Kb] = 2 (g> 2), |*ab) = Ci,o|l, 0} + Ci,q|0, 1). 

For state (p3) the probability of counting k photons is independent of the coefficients Cjv— n,ni 



N- 



(52) 



1. One illustrative example of this kind 



P( k , t) = ^ (^j ^ (-e-* + 4e-^ 2 + <yi - 3) fc 

x exp [- (2x/t) 2 iV 2 (-e~ 7 * + 4e~ 7t/2 +jt-3) 



(53) 



so, all higher moments are determined from the first one, the mean counted photons is a precise measurement since the 
variance is zero and consequently the squared coefficients \CN-n,n\ 2 cannot be determined. For 3> 1 the expression 
for the r-moments (expressed as K, see (48)) of counted photons will be 



~W = N 2r = {K) r . (54) 
(b) Now, let us consider the case when both modes are prepared in the superposition of perfect correlation 



N/2 

|*ab) = ^ c n,n\n,n) . 

n=0 

Again, for qubits, I^ab) = Co,o|0,0) + Ci, i 1 1 , 1) - For the state fl5q) the probability of counting k photons is 



(55) 



P(fc, t ) = - ( ^ ) J2 \C n ,n\ 2 (2nf k (-e-^ + Ae-^ 2 +1 t- 3)* 



2k N/2 



x exp 



n=0 
2 



(2n) 2 (-e-^ + 4e-'' t / 2 +7t-3) 



(56) 



and for jt 3> 1 the moments of counted photons (expressed in terms of k) will be 

N/2 
n=0 

and all coefficients squared moduli can be determined, 

C(2n) = |C„, n | 2 , n = 0,1,2,..., N/2. 

When N — > oo and 

(tanhr)™ 



cosh? 



(57) 



(58) 



(59) 



state ( p5| ) is a two-mode squeezed state (r is the squeezing parameter) used to establish the quantum channel for the 
continuous variable teleportation, as reported in |24(| . So, all terms can be precisely determined since all depend on 
the parameter r that can be inferred from the relations between the coefficients. 

We call attention to the fact that although this procedure permits to determine the coefficient moduli |C m , m | of 
state (EB), it does not allow to distinguish this state from the mixed state 



N/2 



Pab(0) = \C n , n \ 2 \n, n)(n,n\ 
because the same outcomes are inferred, Eq. (p7|), for the r-moments. 



(60) 







V. SUMMARY AND CONCLUDING REMARKS 



We have considered a monitor subsystem M. (an EM mode, C) coupled to a system of interest S, (two interacting 
modes, A and B), where the interacting term in the hamiltonian contains a QND variable of S: the total quanta 
of modes A and B. We proposed a nondeterministic entanglement generation protocol of the two modes, A and B, 
based on the continuous photodetection theory. By counting destructively k photons of the mode C the amount 
of entanglement of the joint state of modes AB, prepared initially in coherent states, can be controlled. Due to the 
dissipative character of the non-ideal photodectection model, nonmaximally entangled states (mixtures) are generated. 
The distribution function for counting photons k and time intervals t is well defined allowing the determination of 
the most probable time for the occurrence for each value of k. 

For the sake of comparison we also investigated the case of an ideal or projective measurement (instantaneous) 
as discussed by von Neumann. The projection of the S — M. state in a photon number operator eigenstate \k) pure 
entangled of A and B, however no control can be done, because the realization of eigenstate cannot be fixed a priori, 
it occurs probabilistically. 

The post-selected counting distribution function allows calculating the moments of the counted photons of C, which, 
for yt ^> 1 are closely related to the moments of the squared number of photons of both modes, N 2 , which is a constant 
of the motion. So, also higher moments of N 2 can be inferred in a nondcmolition measurement. This has a immediate 
use, as probing the state of the AB system by means of the counting distribution function. Examples were given, 
and as expected, the extracted information showed to be useful for partial or total reconstruction of the initial state 
of modes A and B. Since only the squared modulus of the coefficients of the state are given, we remarked that such 
procedure is not able to distinguish systems between pure or mixed states prepared as maximally entangled. 

Although the system here discussed is constituted of EM field modes, the couplings may possibly be realized in 

] or even Bose-Einstein condensates [pfjj. In 
both cases the coupling of the atomic system with a light field (monitor mode) is able to entangle atomic systems. It 
is also possible to probe the atomic system state by photocounting on a light beam. In such cases we expect collisions 
to be important if not restrictive to the method. An obvious extension of the protocol here proposed is to use the 
classical information achieved to control the AB system state in a continuous feedback process |27j . This allows the 
coherent control of state entanglement of systems, and would be, indeed, useful for quantum information processing 



other systems such as vibrational degrees of freedom of trapped ions 1 25 
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figures 

Fig. (1) - M. C. de Oliveira et al. 




yt 



FIG. 1. Photocount probability distribution for initial coherent states for modes A and B, with \a\ 2 = |/3| 2 = 5. Numbers 
above the curves indicate the counted photons, jt is a dimensionless time scale. 



Fig. (2) - M. C. de Oliveira et al. 




FIG. 2. The excess entropy at the short time limit (filled squares) and excess entropy at most probable time of the fc-event 
(crosses). For k = 0, ...10, at yt = 0.0, 0.32, 0.4, 0.46, 0.51, 0.55, 0.59, 0.62, 0.65, 0.68, 0.71. 



Fig. (3) - M. C. de Oliveira et al. 
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FIG. 3. Measure of purity of the fc-event, at same times as in Fig. 2. 
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